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ABSTRACT 

Recent discoveries by the Sloan Digital Sky Survey (SDSS) of four bright z ~ 6 quasars could constrain 
the mechanism by which the supermassive black holes powering these sources are assembled. Here we 
compute the probability that the fluxes of the quasars are strongly amplified by gravitational lensing, and 
therefore the likelihood that the black hole masses are overestimated when they are inferred assuming 
Eddington luminosities. The poorly-constrained shape of the intrinsic quasar luminosity function (LF) 
at redshift ~ 6 results in a large range of possible lensing probabilities. If the LF is either steep, or 
extends to faint magnitudes, the probability for amplification by a factor (i > 10 (and with only one 
image detectable by SDSS) can reach essentially 100%. We show that future observations, in particular, 
either of the current four quasars at the high angular resolution provided by the Hubble Space Telescope, 
or of an increased sample of ~ 20 redshift six quasars at the current angular resolution, should either 
discover several gravitational lenses, or else provide interesting new constraints on the shape of the z ~ 6 
quasar LF. 

Subject headings: cosmology: observations — cosmology: theory — galaxies: high-redshift — 
gravitational lensing — quasars: general - black hole physics 



1. INTRODUCTION 

The Sloan Digital Sky Survey (SDSS) has recently dis- 
covered several extremely bright and distant quasars with 
redshifts z ~ 6, whose luminosities are a result of gas ac- 
cretion onto supermassive black holes (BHs). These rare 
objects have luminosities at the far bright tail of the quasar 
luminosity function (LF) , corresponding to ~ 5 a peaks in 
the density field in models in which quasars populate dark 
matter halos (Haiman & Loeb 2001; Fan et al. 2001b). By 
assuming that the observed luminosities equal the limiting 
Eddington luminosities of the BHs, the masses of all these 
BHs are inferred to be a few x 10 9 M Q . The mere existence 
of such massive BHs at so early a stage in the evolution 
of the universe should provide insight into the formation 
and growth of supermassive black holes (Haiman & Loeb 
2001). 

The above estimates for the masses of the BHs asso- 
ciated with the SDSS quasars scale directly with the in- 
ferred luminosities, and thus with the observed fluxes of 
the quasars. However, the apparent fluxes can be increased 
by gravitational lensing. If a quasar is gravitationally 
lensed and this is not taken into account in the mass es- 
timate, the quasar's intrinsic luminosity, and therefore its 
BH mass, will be overestimated. Previous considerations 
of the probability that these SDSS quasars are lensed in- 
dicate that while the a priori lensing probability is small, 
the magnification bias is poorly known and can be quite 
large (Wyithe & Loeb 2002). If massive BHs grow from 
stellar-mass seeds by accreting gas at the Eddington rate, 
they reach masses of ~ 10 9 M Q in ~ 10 9 years, a time 
span comparable to the age of the universe at z ~ 6 in the 
current "best-fit" cosmologies. 
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Since the presence of such massive BHs can only be 
marginally accommodated into structure formation mod- 
els (Haiman & Loeb 2001), it is interesting to ask whether 
all of the SDSS quasars may be strongly lensed. In this 
paper, we compute the expected probabilities for lensing 
of z ~ 6 quasars by intervening galaxies. We quantify 
how the probability depends on the shape of the assumed 
quasar LF, and in particular, we focus on the question: 
can the expected probabilities reach values near unity? We 
show that current constraints on the z ~ 6 quasar LF are 
sufficiently weak that the probabilities for strong lensing 
can indeed be close to 100%. We will utilize this result to 
"invert" the problem, and demonstrate how the intrinsic 
LF can be constrained using current and future detections 
(or absence) of lensing events with multiple images. 

The rest of this paper is organized as follows. In § 2, 
we describe our model for the population of gravitational 
lenses. In § 3, we derive the resulting intrinsic "a priori" 
lensing probabilities along a random line of sight to red- 
shift z ~ 6. In § 4, we discuss the effects of magnification 
bias, and the expected a posteriori lensing probabilities for 
the z ~ 6 quasars in models with different LFs. In § 5, 
we quantify constraints on the quasar LF that can be de- 
rived from current and future searches for lensing events. 
Finally, in § 6, we discuss our main results and then sum- 
marize the implications of this work. Throughout this pa- 
per, we adopt a flat cosmological model dominated by cold 
dark matter (CDM) and a cosmological constant (A), with 
fi m = 0.25, tt b = 0.04, and Cl A = 0.75, a Hubble constant 
Hq = 70 km s _1 , an rms mass fluctuation within a sphere 
of radius 8 h^ 1 Mpc of er 8 = 0.9, and power-law index 
n = 1 for the power spectrum of density fluctuations. We 
also adopt the cosmological transfer function from Eisen- 
stein & Hu (1999). 
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2. MODEL DESCRIPTION 

The a priori probability of gravitational lensing depends 
primarily on the comoving number density of lenses and 
the lensing cross section of each lens. Present theories 
agree that galaxies are created by the condensation of gas 
in dark matter halo cores, and the abundance of the most 
massive halos depends on the overall amplitude of mass 
fluctuations. We use the halo mass function of Jenkins et 
al. (2001), who measure dn/dM, the number of halos per 
unit comoving volume and per unit mass, from TV-body 
simulations. 

We employ a combination of two separate lens density 
profiles. The simpler of the two halo density profiles we 
consider is the singular isothermal sphere (SIS) density 
profile. The SIS profile is characterized by its depen- 
dence on the one-dimensional velocity dispersion, which 
is supported by observations of flat rotation curves, and 
its density profile is p(r) = (T^/2-KGr 2 , where a v is the 
velocity dispersion. The steep inner density slope of the 
SIS profile, however, renders it a poor fit to the more 
massive, and more extended halos. As a result, Navarro, 
Frenk, & White (1997; hereafter NFW) proposed a "uni- 
versal" mass density profile for dark matter halos found 
in numerical simulations. The NFW density profile is 
p(r) = p s r 3 /r(r + r s ) 2 , where p s and r s are constants. 
Although it matches the SIS model at intermediate radii, 
it branches off for large and small radii. 

Because of the dynamics of cooling gas in halos (Keeton 
1998; Porciani & Madau 2000; Kochanek & White 2001; 
Li & Ostriker 2002), we use a combination of the SIS and 
NFW profiles to describe the halos that contribute to the 
quasar lensing probability. Gas falling into very massive 
halos is shock-heated to high temperatures (T > 10 7 K) 
and cannot cool within a Hubble time. The gas therefore 
remains in hot, extended gaseous halos around groups or 
clusters of galaxies which have shallow inner density pro- 
files. For less massive halos, in contrast, the virial temper- 
ature is low enough that the gas can cool and collapse into 
galaxies with steep inner density slopes. This difference in 
density slopes implies that a single density profile cannot 
describe all lensing objects in the universe. 



As a solution, we use the SIS profile [dp/dr 
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less massive halos which have steep inner density slopes 



and the NFW profile (dp/dr 



for r <C r s ) for more 



massive halos which have shallow inner density slopes. We 
set the transition between these two profiles at a dark mat- 
ter halo mass of 10 13 M©, the mass at which typically half 
of the dark halos have cooled the majority of their gas 
content, as given by the halo cooling probability distribu- 
tion in Kochanek & White (2001). While this value of 
the cutoff has a physical basis, perhaps more importantly 
the predicted splitting angle distribution for lensing events 
with multiple images matches observations (Keeton 1998; 
Kochanek & White 2001; Li & Ostriker 2002). Therefore, 
in our calculations, we used the SIS profile for halo masses 
less than 10 13 M© and the NFW profile for halo masses 
greater than 10 13 M©. 

NFW halos are very inefficient lenses, with a cross- 
section per unit mass > 100 times smaller than for SIS 
halos. Although the high mass halos are responsible for 
most lensing events with large (> 5") image splittings, 
these comprise only a small (<~ 1%) fraction of all lensing 



events of interest. As a result, in practice, we find that 
we can safely ignore NFW halos, and use a model where 
halos with masses below 10 13 M Q are assumed to have SIS 
profiles, and halos with masses above 10 13 M© are simply 
excluded from the calculations. Our results for the in- 
trinsic lensing probability are relatively insensitive to the 
choice of this mass cutoff: a factor of two higher/lower 
value would increase/decrease the intrinsic lensing proba- 
bility to z ~ 6 by 20%. The effect on the a posteriori 
lensing probability is similar in models where the proba- 
bility is relatively low (< 50%), but the mass cutoff has a 
smaller effect on the a posteriori lensing probability if this 
probability is high (see discussion below). 




Fig. 1. — The a priori probability as a function of lens redshift 
that the z = 5.99 SDSS quasar was magnified by at least a factor 
Mmin = 2, but has not produced two separately detectable images 
for SDSS. The middle curve in the figure shows the probability for 
multiple images with an image splitting of A9 < 1". The bottom 
curve shows the probability for image splittings of A0 > 1", but 
subject to the constraint that the fainter image must be below the 
SDSS detection limit (this corresponds to 2 < fi < 6.172 for this 
source). The top curve shows the total lensing probability. 



3. A PRIORI LENSING PROBABILITY 

The a priori probability that a quasar is lensed by an 
intervening galaxy is given by the product of the number 
of lenses per unit volume, (1 + z) 3 dn/dM, the galaxy 
lensing cross section, , and the path length of light from 
the quasar, c (dt/dz) dz, where dt/dz = —1/(1 + z)H(z), 
integrated over all possible lens masses Ml and redshifts 
Zl- Therefore, the a priori probability that a source at 
redshift z s is lensed and magnified by a factor > p is 

p(p,z s ) = / c L o- L {> p;M L ,z Ll z s ) 

Jo Jo tiyz L ) 

dn 

x^rj(M L ,z L )dz L dM L , (1) 
dM 

where the Hubble parameter at redshift z in a flat uni- 
verse is H(z) = H [n m (l + zf + fi A ] 1/2 , dn/dM is the 
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Fig. 2. — The a priori probability for the z = 5.99 SDSS quasar 
to be magnified by at least a factor of two, as a function of the 
maximum allowed splitting angle AO. The dashed curve includes 
all halos, and assumes that each halo is described by an SIS profile. 
The solid curve excludes halos above 1O 13 M0. Including halos above 
this mass, described by NFW profiles, would yield a curve that is 
nearly identical to the solid curve. For splitting angles of AO > 1" 
we require that the fainter image must be below the SDSS detection 
threshold. This excludes lensing events with fi > 6.172, and causes 
a ~ 10% decrease in the lensing probability for A8 > 1". 



comoving halo mass function from Jenkins (2001), and <jl 
is the cumulative lensing cross-section in the lens plane for 
lensing amplification by a factor > jj, by a SIS (see Schnei- 
der, Ehlers & Falco 1999). Note in particular that the 
cross section scales with the magnification as <tl oc 
for fi > 2, and as ctj, oc (/i — 1)~ 2 for <2. 

Fig. 1 shows the a priori probability dp(fi = 2, z s = 
5.99) /dzL as a function of lens redshift that the z = 5.99 
quasar SDSSp J130608. 26+035636.3 (magnitude z* = 
19.47) was magnified by at least a factor /j m j n = 2, but has 
not produced two images separately detectable by SDSS. 
The middle curve in the figure shows the lensing probabil- 
ity for multiple images with an image splitting of A9 < 1", 
approximately the angular resolution of the SDSS (this 
corresponds to a redshift-dependent upper-mass cutoff of 
~ 10 12 M Q for the lensing halos). The bottom curve shows 
the probability for image splittings of A9 > 1", but sub- 
ject to the constraint that the fainter image must be be- 
low the SDSS detection limit of z* = 20.2. Note that 
for an SIS lens, multiple (i.e., two) images are formed for 
[i > 2, with a flux ratio ^1/^2 = (Mtot + 2)/(/i t ot - 2). 
Thus, the flux ratio is close to unity for high magnifica- 
tions, and increases monotonically for lower total magnifi- 
cations. In practice, the constraint that only one image be 
detectable simply imposes a maximum total magnification 
of Mmax = 2(10°- 4Am + l)/(10°- 4Am - 1). This source is 
Am = 0.73 mag above the detection threshold and thus 
Mmax ~ 6.172. The top curve in Fig. 1 shows the total 



lensing probability. As the figure shows, the total prob- 
ability is very small, of order 0.1 percent, and reaches a 
maximum near a lens redshift of zl = 1.5. 

Fig. 2 illustrates how the total lensing probability, p(p, > 
2,z s — 5.99), depends on the maximum allowed image 
separation. For small splitting angles the probability in- 
creases steeply, whereas for large angular separations it 
flattens out. The dashed curve is for all halos, while the 
solid curve only includes halos less massive than 10 13 M Q . 
As discussed above, for the purposes of gravitational lens- 
ing, halos more massive than 10 13 M Q are better described 
by NFW profiles, which have a much smaller cross-section 
for lensing than halos with SIS profiles. Including the 
M > 10 13 M Q halos using NFW profiles would result 
in a curve that is nearly identical to the solid curve. In 
both curves shown in this figure, for splitting angles of 
A9 > 1", we require that the fainter image must be below 
the SDSS detection threshold. As discussed above, this ex- 
cludes lensing events with /x > 6.172, and causes a ~ 10% 
decrease in the lensing probability for A9 > 1". Note that 
according to Figure 2, ~ 75% of all lensing events have 
image separations below 1", but only ~ 10% have image 
separations below 0.1". 

Our results in this section agree to within a factor of 
two with previous estimates of the lensing probability and 
its dependence on lens redshift, e.g. by Kochanek (1998) 
and by Barkana & Loeb (2000). We conclude that the 
a priori lensing probability for az~6 SDSS quasar along 
a random line of sight is small, and dominated by lenses 
between redshifts 1 < zl < 2. 

4. A POSTERIORI LENSING PROBABILITY 

The a priori lensing probability calculated in §3 reflects 
the probability that lensing causes magnification by a fac- 
tor of at least (i along a random line of sight. However, we 
are interested in the a posteriori probability that a given 
SDSS quasar has been lensed. Since the SDSS can only 
detect z ~ 6 quasars that are at the bright-end tail of the 
LF (if extrapolated from lower redshift), the magnification 
bias, defined as the ratio between the a posteriori and the 
a priori probability, could well be significant. 

We determine the a posteriori lensing probability from 
the observed quasar LF. If we define <J>i nt (L) as the intrin- 
sic space density (assuming no lensing) of quasars with 
luminosity L ± dL/2, then, given the intrinsic differen- 
tial lensing probability dp/d^i, the space density of quasars 
that will be observed at luminosity L b s ±dL b s /2 and red- 
shift z s , and that are magnified by at least a factor ^ m j n 
is 



< I , obs(iobs 1 /^min ) 



dfl^$ int (L ohs / fj,)^-. (2) 



Similarly, the total space density of quasars observed at 
the same luminosity and redshift, regardless of their ampli- 
fications, is $ bs(£ bs, 0) and the a posteriori probability 
that a quasar with observed luminosity L \, s is lensed by 
a minimum factor /x m i„ is 

$obs(£obs 



(3) 



$obs(£obs,0) 

The magnification bias and the a posteriori lensing prob- 
ability depend strongly on the shape of the quasar LF. If 
the LF is steep or extends to faint luminosities, the mag- 
nification bias can be strong. Note that the magnification 
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bias depends only on the shape ol the LF, not on its nor- 
malization. For the purpose of calculating P(> Mmin), we 
need, in principle, to extend the intrinsic lensing proba- 
bility distribution to cover the full range of amplifications, 
including de-amplification, < ji < 2. We shall use the 
following form of dp/dji: 



dp 
dji 



8t 2 /m 3 (P > 2) 

T 2 /(fi - l) 3 (/i < n < 2) 

r 2 /(Mo - I) 3 (1 < It < Mo) 

(/x < 1) 



(4) 



Here r 2 = / 2 °° d/j,(dp j 'd[i) is the total optical depth for 
forming double images (which depends on the source red- 
shift, and is obtained from equation 1). The top two rows 
in equation (4) describe an SIS lens as usual, with a nor- 
malization that ensures a smooth transition from the case 
when the source is doubly imaged (but unresolved) to the 
singly imaged case. In order to correctly model the prob- 
ability distribution near fi = 1, we would have to allow 
for magnifications smaller than unity, which would require 
modeling the large-scale density distribution of the uni- 
verse. Instead, we have imposed a critical magnification 
jUo ~ 1, below which the probability is set to a constant. 
The value of /zq is chosen to be 



Mo 



1 



3r 2 



(5) 



1. Note that 



2-t 2 

which enforces the normalization J Q C 

this simple characterization does not strictly conserve flux, 
since the mean magnification J °°rf/x^/i > 1. However, 
in practice, we find that (Xq and the mean magnification 
are close to unity, and our results are insensitive to our 
treatment of dp/dfj, near fj, = 1. 

To calculate the a posteriori lensing probability of the 
high rcdshift SDSS quasars, we also need to know their in- 
trinsic LF. We here assume that the intrinsic z ~ 6 quasar 
LF can be described by the double power law form given 
in Pei (1995): 

$int(i) = (L/L^fil/L^ - (6) 

The same form was used by Madau, Haardt, & Rees (1999) 
and more recently by Wyithe & Loeb (2002) , with different 
sets of parameters. The LF depends on four parameters: 
the normalization <!>», the faint-end slope fii, the bright- 
end slope fih, and the characteristic luminosity at which 
the LF steepens. The a posteriori lensing probability is 
sensitive to the last two parameters, (3h and L* (and is 
strictly independent of the normalization <&*). 

Through different combinations of slopes and breaks, 
we present seven possible LFs at redshift z = 6, which are 
illustrated in Fig. 3. All models have = 1.64 (Pei 1995) 
and the normalization $* has been adjusted so that each 
curve passes through the Fan et al. (2001b) data point 
corresponding to the space density implied by the four 
z ~ 6 quasars. The curves in the upper panel all have 
fih ~ 3.5, but varies from « 10 10 Lb.q for the Pei 
(1995) model (solid curve) to w 10 12 Lb.q for the model 
of Madau et al. (1999). The curves in the bottom panel 
all have L* 10 115 Lb,q, but have different values of fih- 
The parameters of the LFs are summarized in Table 1. 



Table 1 

Parameters of the z — 6 quasar luminosity functions. 
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Fig. 3. — The seven different LFs used to compute lensing prob- 
abilities. The LFs are all described by a double power law with 
the same faint— end slope of /3; = 1.64, and are shown at redshift 
z = 6. The data point is from Fan et al. (2001) and represents 
the space density of luminous quasars at z Kb 6, inferred from the 
discovery of the four SDSS quasars (log(L<E>/Mpc ) = —10.5 at 
log(Ls/LB,Q) = 13.3). The upper panel demonstrates the effect of 
varying the break of the LF. From top to bottom the four LFs have 
log(L*/LB,0) ~ 10.1, 10.9, 11.5, and 11.9 respectively. Apart from 
the normalization <!>*, which has been adjusted to ensure consistency 
with the SDSS data point, the curves correspond to the model of 
Pei (1995) (solid), the Pei (1995) model with a 6 times larger L* 
(dotted), the model of Wyithe & Loeb (2002) (short dashed), and 
Madau ct al. (1999) model (long dashed). The Wyithe & Loeb 
(2002) model (short dashed) has (3^ = 3.43, while the other mod- 
els have f3h = 3.52. The lower panel demonstrates the effect of 
varying the bright-end slope f3 h from 2.5 to 4.5. All three models 
have log(L*/LB,0) ~ 11.5. See Table 1 for a summary of the LF 
parameters. 



The values are given at redshift z = 6; we allow for slightly 
different values for the characteristic luminosity over the 
range 5.8 < z < 6.28 by including the redshift dependence 
of L* as parameterized in each case by Pei (1995), Madau 
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ct al. (1999), and Wyithe & Loeb (2002). 

The lensing probability depends on the intrinsic LF, 
which is not directly observable. In using the Fan et al. 
(2001b) data point as a constraint on the intrinsic LF, we 
thus implicitly assumed that lensing does not have a large 
overall effect on the LF. If magnification bias is significant, 
then this assumption breaks down for the bright-end of the 
LF. Thus, to match the observations, the normalizations of 
the intrinsic LFs should in fact be lowered by a factor equal 
to the magnification bias for the observed quasars. How- 
ever, for the purpose of computing the a posteriori lensing 
probability this is irrelevant, since the magnification bias 
is independent of the normalization of the intrinsic LF. 
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Fig. 4. — A comparison of the intrinsic and observed LFs for 
our most extreme model, which has an intrinsic bright-end slope of 
/3ft = 4.5 (see dashed curve in lower panel of Fig. 3). The lensing 
probability at z = 6 in this case is high, and the apparent LF is 
increased by over an order of magnitude at the luminosity of the 
SDSS quasars. The apparent slope is flattened to fi as 3, and the 
model is therefore consistent with the upper limit on slope of the 
apparent LF from the SDSS ((3 < 4.3 at 99 percent confidence, Fan 
ct al. 2001). 



The SDSS observations provide an upper limit on the 
slope of the bright-end of the observed LF: Ph,obs < 4.3 at 
99 percent confidence (Fan et al. 2001a, b). If magnifica- 
tion bias is important, then this constraint will always be 
satisfied since lensing by an isothermal sphere will result in 
a bright-end slope of [3h,obs = 3.0 (the intrinsic differential 
lensing probability dp/d/j, oc /z -3 for /i > 2). This is illus- 
trated in Fig. 4, which shows the intrinsic and observed 
LFs for our most extreme model, which has (3h = 4.5 
(dashed curve in the lower panel of Fig. 3). Since the 
magnification bias increases with (3h, observations cannot 
directly constrain the steepness of the intrinsic, bright-end 
slope of the LF. However, we will see below that the slope 
can be constrained indirectly, by measuring the fraction of 
multiply imaged quasars. 

Fig. 5 shows the a posteriori probabilities for the z = 




ooi 1 1 1 1 1111 1 1 1 1 1 1 1 

1 10 100 



Fig. 5. — Lensing probability as a function of minimum magnifi- 
cation fi for the z = 5.99 SDSS quasar and the seven different LFs 
shown in Fig. 3. Note that the total lensing probability is sensi- 
tive to variations of both the slope 0^ and break L*, and can be 
high. For fi > 6.172, both lensed images would appear brighter than 
SDSS's threshold, and we impose the condition A6 < 1" on their 
splitting angle to avoid detection by SDSS. On the other hand, the 
fainter image is too faint to be detectable, and we allow lensing with 
arbitrarily large splitting angles for fi < 6.172. 



5.99 SDSS quasar to be magnified by at least a factor /j,, 
subject to the constraint that only one image is detectable 
by SDSS. The small discontinuity in the probabilities at 
\i = 6.172 occurs because for \i < 6.172 the fainter im- 
age is too faint to be detectable, and we do not impose 
any constraint on the splitting angle. On the other hand, 
for \i > 6.172 both images would be brighter than SDSS's 
threshold, and we therefore impose the condition A9 < 1" 
on the splitting angle. The seven curves shown in Fig. 5 
correspond to the seven intrinsic LFs in Fig. 3. From the 
figure it is clear that a wide range of probabilities is pos- 
sible for a given minimum magnification. In particular, 
because of the lack of constraints on the intrinsic LF, they 
can be significantly higher than the maximum value of 
- 30% found by Wyithe & Loeb (2002), who considered a 
more restricted range of LFs (the short dashed curve in the 
upper panel in Fig. 3 represents their most extreme model 
that results in their maximum lensing probabilities). For 
example, if L* is small, as for the Pei (1995) model (solid 
curve in upper panel), then the probability that the quasar 
has been magnified by more than a factor of ten is about 
sixty percent. The lensing probability also becomes very 
large if the bright-end slope is steep: P{[i > 10) is close 
to unity if [3h = 4.5 (dashed curve in the lower panel). In 
Table 2, in the last line, we list the lensing probability for 
each of the four quasars in this model. 
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Table 2 

Parameters and lensing probabilities of the four z « 6 SDSS quasars 8. 
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P((U > Mm,, 


) [AO > 1' 


',/?)» 


= 4.5] 


0.220 
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P(p > 10) 


[A6» < 1' 


,A, 


= 4.5] 


0.877 


0.939 


0.862 
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a The first two lines list the redshift and magnitude of each source. The 3 rd and 4 th 
lines give the optical depth to lensing with amplification ^ > 2 with and without including 
events whose splitting angles are larger than l". The 5 th line lists the minimum total 
magnification for which both images arc above the SDSS magnitude threshold of z* — 20.2. 
The last two lines list the lensing probabilities, with magnification bias taken into account, 
in our most extreme model with a bright-end quasar LF slope of 0h — 4.5. The 6 line 
gives the probability for a event that produces two separately detectable images for SDSS. 
The last line gives the probability for an unresolved lensing event with total magnification 
ft > 10. 



5. CONSTRAINING THE SHAPE OF THE LUMINOSITY 
FUNCTION 

In the previous section we have demonstrated that it is 
possible that most of the z ~ 6 SDSS quasars are strongly 
magnified by gravitational lensing, and that the expected 
lensing probability depends very strongly on the shape of 
the intrinsic quasar LF. Because of the unknown magnifi- 
cation bias, surveys that probe only the bright-end of the 
LF (L > L*), cannot directly measure the intrinsic LF. 
However, it is possible to constrain the intrinsic LF, in 
particular the characteristic luminosity L* and the bright- 
end slope 0h, by measuring the fraction of multiply imaged 
quasars down to a fixed limiting angular resolution. 

The SDSS survey has so far discovered four z ~ 6 
quasars (Fan et al. 2001), all of which appear to be sin- 
gle images at the angular resolution (<~ 1") and detection 
limit (z* — 20.2) of the survey. The parameters of these 
quasars are listed in Table 2. Lensing by an isothermal 
sphere always yields a double image if the total magnifica- 
tion is greater than two. Thus, if these object are strongly 
lensed, then the extra image must have escaped attention 
either because the splitting angle was too small for SDSS 
to resolve the images or because one of the images was 
too faint to be detected. In the previous section, when 
we computed the a posteriori probabilities for lensing, we 
took these limitations of the survey into account. 

In this section, we relax these limitations, and compute 
the complementary quantity: the probability of lensing 
events that would produce multiple images, both of which 
would have been detectable by SDSS (i.e., we require that 
both images be above the detection threshold, and sepa- 
rated by > 1"). The lack of such a detection can be used 
to constrain models for the LF. In particular, for a given 
quasar LF, we can compute the likelihood that at least one 
of the four SDSS quasars would exhibit two images. This 
is given by P = 1 - (1 -Pi) x (1 -P 2 ) x (1 -P 3 ) x (1 -P 4 ), 
where Pi is the probability for the i th quasar to be doubly 
imaged. Wc consider variations in the two parameters 
and 0h of the quasar LF, and compute P as a function of 
and (3h- As an example, in Table 2, in the 6 th line, we 
list the probability for each quasar to be doubly imaged in 
our most extreme model for the quasar LF. 



5 - 



5 - 
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Fig. 6. — Constraints on the parameters of the redshift z = 
6 quasar LF from the non-detection of gravitational lenses. The 
solid curves show likelihood contours for ruling out a given LF, if 
no gravitational lenses are found. The upper left panel shows the 
la constraint from the existing four SDSS quasars, assuming that 
double images arc resolvable if they are separated by > 1", and are 
both above the magnitude threshold of z* = 20.2. The constraint is 
weak. The other three panels demonstrate how the constraint can 
be improved if no lenses are found despite (i) increasing the angular 
to 0.1" (upper right panel); (ii) increasing the sample size to a total 
of 20 quasars (lower left panel); or (iii) both (lower right panel). 
The 3cr contour is off the scale of the plots, except in the lower right 
panel; the 2<j contour is off the scale in the upper left panel. The 
dashed curves show contours of fixed ionizing background T_i2 = 
0.1 and T_i2 = 2. 



In the upper left panel of Fig. 6, we show the constraint 
on the intrinsic LF resulting from the non-detection of any 
multiple images in the four quasars reported in Fan et al. 
(2001). The solid curve corresponds to the constant P = 
0.68, i.e., models that lie above the solid curve are excluded 
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at the 1 a level. The analogous 2 and 3 a contours are off 
the scale of the plot. The SDSS observations do not yet 
place strong constraints on the intrinsic LF. This can be 
understood as follows: the intrinsic probability for lensing 
with multiple images to z = 6 when all lensing halos are 
considered (up to a mass 10 13 M©, see § 3) is t 2 = 0.00326, 
but when splittings of > 1" are excluded, the probability 
is r 2 = 0.00247. In other words, only ~ 25% of all lensing 
events would be resolved at > 1". Even for an arbitrarily 
large magnification bias, SDSS would only have a ~ 25% 
chance of resolving the multiple images of a given quasar 
(see Tabic 2). Note that the z = 6.28 source is very close 
to the SDSS magnitude threshold, and the fainter image 
would escape detection unless the total magnification is 
jU > 29. As a result, the probability for producing an 
additional detectable image for this source is smaller than 
for the other three quasars. 



aa. 



4 ^ 
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Fig. 7. — Likelihood contours as in Fig. 6, but for a single z = 6, 
z* = 20 quasar. The curves in each of the four panels show con- 
stant probabilities for a lensing event that is detectable at different 
angular resolutions (the upper right panel assumes infinite resolu- 
tion, i.e., it includes all multiple lensing events). Each panel shows 8 
curves, corresponding to the following probabilities (bottom to top, 
in percent): (i) Upper left panel: 1, 2, 4, 10, 20, 23, 24, 24.18; (i) Up- 
per right panel: 4, 8, 15, 40, 80, 95, 99.3, 99.9; (i) Lower left panel: 
2, 4, 9, 20, 35, 43, 44.3, 44.55; (i) Lower right panel: 3, 6, 13, 33, 60, 
72.5, 75.4, 75.86. Note the "saturation" of the probabilities near the 
top contour in each panel, caused by the finite fraction of lensing 
events that produce splittings above a given angular resolution. 



In the remaining three panels of Fig. 7 we illustrate how 
future observations would constrain the LF if no multiple 
images were detected. The bottom left panel shows the 
results for 20 objects with properties identical 2 to the four 
already discovered quasars. If SDSS were to discover 20 
z ~ 6 quasars, and none of them showed a second image, 

2 In increasing the sample size, we simply replicated the current sam- 
ple of quasars five times. Since the lensing probabilities are similar 
for quasars near the detection threshold and near z = 6, our results 
are not sensitive to this particular choice. 



then the more extreme intrinsic LFs plotted in Fig. 3 (the 
Pei 1995 and the j3u — 4.5 models) would be ruled out at 
the 2 a level. 

The upper right panel shows that nearly identical con- 
straints can be obtained by observing the four already 
known quasars with an angular resolution of 0.1", a reso- 
lution typical for observations with the Hubble Space Tele- 
scope (for the cases of a maximum splitting angle of 0.1" 
we assumed that sufficient depth is reached that both im- 
ages arc detectable, and did not impose a magnitude cut, 
or a constraint on the flux ratio) . Finally, the bottom right 
panel shows that observations of 20 quasars at 0.1" resolu- 
tion would enable us to put very strong constraints on the 
intrinsic LF, with the current "ball-park" LFs, obtained 
by extrapolations from lower redshifts, ruled out at the 
~ 3 a level. 

In these examples we assumed that no multiple images 
are detected. Similar calculations can of course be done 
if future observations would detect one or more multi- 
ple images. In that case the likelihood contours can be 
closed and the parameters of the LF would be constrained 
both from above and from below. To enable interpreta- 
tion of any outcome of a future search for lenses, Figure 7 
shows probability contours for multiple lensing for a sin- 
gle z — 6, z* = 20 quasar. The curves in each of the 
four panels show constant probabilities for a lensing event 
detectable at four different angular resolutions (the upper 
right panel assumes infinite resolution, i.e., it includes all 
multiple lensing events). Each panel shows 8 curves, corre- 
sponding to the following probabilities (bottom to top, in 
percent): (i) Upper left panel: 1, 2, 4, 10, 20, 23, 24, 24.18; 
(i) Upper right panel: 4, 8, 15, 40, 80, 95, 99.3, 99.9; (i) 
Lower left panel: 2, 4, 9, 20, 35, 43, 44.3, 44.55; (i) Lower 
right panel: 3, 6, 13, 33, 60, 72.5, 75.4, 75.86. Note the 
"saturation" of the probabilities near the top contour in 
each panel, caused by the finite fraction of lensing events 
that produce splittings above a given angular resolution. 

Varying the intrinsic LF does not only affect the lensing 
probability, but also changes the total amount of emit- 
ted quasar light, including the ionizing emissivity. Since 
the ionizing background can in principle be measured, 
this can provide an additional constraint on the quasar 
LF. Observations of the mean absorption in the Lyman- 
a forest have been used recently to infer the ionizing 
background at z ~ 6 (Becker et al. 2001; Cen & Mc- 
Donald 2001; McDonald & Miralda-Escude 2001; Fan ct 
al. 2002), yielding ionization rates per hydrogen atom of 
r_ 12 = T/lO-^s- 1 « 0.1. The dashed curves in Fig. 6 
show contours of fixed quasar emissivity, corresponding to 
r_i2 = 0.1 and 2 respectively. Our method for computing 
the ionization rate from the LF is described in more detail 
in the Appendix below. At face value, LFs that lie above 
the lower dashed curve over-produce the ionizing back- 
ground; those above the upper dashed curve by a factor of 
> 20. However, in practice this result has large uncertain- 
ties. First, although the adopted quasar LF yields the to- 
tal ionizing emissivity, the prediction of the ionizing back- 
ground involves the mean free path of ionizing photons, 
which is highly uncertain at the redshifts of interest here. 
A small mean free path would reduce the background. Sec- 
ond, a smaller fraction of ionizing photons might escape 
from the dense surroundings of high-redshift quasars than 
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from their lower-redshift counterparts. 

6. DISCUSSION AND CONCLUSIONS 

In this paper, we computed the expected lensing prob- 
abilities of z ~ 6 quasars by intervening galaxies. The 
a posteriori probability that a z ~ 6 quasar is strongly 
lensed (by a factor /i > 10 in flux enhancement), but with- 
out producing two images detectable in the Sloan Survey, 
can be very large, depending on the shape of the LF. The 
LFs we use are consistent with other available observa- 
tional constraints, yet these constraints still allow for a 
wide range of possibilities. 

Though an observed bright-end slope of ph — 4.3 could 
be ruled out at the 99 percent confidence level (Fan et 
al. (2001)), lensing alters the slope so that the apparent 
LF will have a slope close to ~ 3.0 if magnification bias 
is important. As illustrated in Fig. 4, even the steepest 
intrinsic LF we consider, with a bright end slope of 4.5, is 
consistent with the SDSS limit on slope. This will be true 
for arbitrarily steep slopes, provided that the characteristic 
luminosity is not increased. 

In a similar probability calculation, Wyithe & Loeb 
(2002) concluded that the observed probability reaches 
30% that a single z ~ 6 quasar is lensed and magnified 
by a factor of 10 or more. If we restrict our analysis to 
the choices of LFs considered in their paper, we agree with 
this result. However, we here use a wider range of quasar 
LFs and conclude that the lensing probability can reach 
essentially 100%. As a result, the current observations are 
consistent with all four z ~ 6 SDSS quasars being strongly 
lensed. 

If the probability that all four quasars are lensed were 
high, then this would alleviate the problematic time con- 
straints on assembling supermassive black holes at the ear- 
liest stages in the evolution of the universe . However, in 
a separate analysis, based on the apparently large size of 
the ionized region around the SDSS quasar at z = 6.28, 
Haiman & Cen (2002) place a strong constraint on the 
lensing magnification of this source, and find /i < 5. This 
result depends primarily on the assumption that the source 
is embedded in a neutral (rather than ionized) intergalac- 
tic medium. 

We have also considered the probabilities for lensing 
events that would have been detectable by SDSS. The 
lack of such detections can be used to place constraints 
on the quasar LF. Although constraints from the current 
four quasars are mild, the situation is likely to improve. 
As illustrated in Fig. 6, increasing the SDSS sample from 
4 to 20 objects, something that is likely to happen over 
the next few years, would allow one to rule out interesting 
quasar LF models. The expected probability of detecting 
multiple images is also sensitive to the angular resolution 
of the observations. In terms of the constraints on the 
shape of the intrinsic LF, observing four quasars at a reso- 
lution of 0.1" is roughly equivalent to observing 20 objects 
at the resolution of the SDSS (1"). 

Hence, upcoming observations, in particular with the 
Hubble Space Telescope, should reveal whether a significant 
fraction of the z ~ 6 quasars is lensed, and will allow us 
to place strong constraints on their intrinsic LF. 
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APPENDIX 

Here we discuss constraints on the high-redshift quasar 
LF from their integrated ionizing radiation. Given an 
LF, we compute the total comoving emissivity as ve u = 
J L$(L)dL. This gives ve v at rest-frame 4400 A; how- 
ever, with a spectral slope of e v oc t/ _1 this is independent 
of wavelength, and has the same value at 912 A. The back- 
ground ionization rate per hydrogen atom T is then given 

by 



no- H {l + zfve v _j 

1 — — A m fp ; S 



4ir 



(7) 



where we have assumed e v oc v~ x , er# = 6.3 x 10~ 18 cm~ 2 
is the hydrogen ionization cross section at the Lyman limit 
(and we assumed o~h{v) oc v 3 above the threshold), hp is 
the Planck constant, i^j = 3.29 x 10 15 Hz is the frequency 
of a Lyman limit photon, and A m f p is the mean free path 
of ionizing photons. We also define r_i 2 = T/10~ 12 s _1 . 
We assume that the effective mean free path corresponds 
to a redshift interval of Az = 0.17 at z = 3 (see Haardt & 
Madau 1996 and Steidel, Pettini & Adelberger 2001), and 
scales with redshift as A m f p oc (1 + z)~ 6 towards higher 
redshift (Cen & McDonald 2002). Based on recent mea- 
surements (e.g. Cen & McDonald 2002, see also Becker et 
al. 2001; McDonald & Miralda-Escude 2001 and Fan et al. 
2002), we require r_i 2 < 0.1 s _1 at z=6. We find that in 
some of our models, at z = 6, this bound is violated by up 
to a factor of 13: r_ 12 = 1.3 (Pei LF); r_ 12 = 0.39 (WL 
LF with a bright-end slope of j3 h = 4.5); r_ 12 = 0.08 (Pei 
LF with 6 times higher L»). More generally, contours of 
constant r_i 2 = 0.1 and T_i 2 = 2 are shown in Figure 6. 

The apparently large ionizing emissivities imply that 
in our LF with the lowest L», the background is over- 
produced, although this result is subject to several uncer- 
tainties : (1) the escape fraction of ionizing photons may 
be only ~ 10% at z — 6; (2) the mean free path evolves 
more steeply than (1 + z)~ 6 ; (3) the faint-end slope is 
considerably shallower than fli = 1.64, or (4) there can be 
significant systematic errors on the observational determi- 
nation of r_i2, since it relies on cosmological simulations 
for modeling the temperature and density field. 
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